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ABSTRACT 


Convergence  zone  (CZ)  prediction  models  are  developed 
based  on  acoustic  ray  tracing  theory  as  applied  to  linearly 
segmented  sound  velocity  profiles  (SVP) .   The  models  were 
developed  into  three  calculator  programs,  two  for  CZ  range 

predictions  under  different  source  and  receiver  depth  con- 
ditions and  one  for  CZ  gain  and  transmission  loss  (TL)  pre- 
dictions.  The  performance  of  the  models  as  programmed  on 
Hewlitt-Packard  HP-67  or  HP-97  programmable  calculators  was 
compared  to  the  Fast  Asymptotic  Coherent  Transmission  (FACT) 
model  which  is  based  on  similar  but  more  elaborate  theory 
and  which  is  designed  for  use  on  large  digital  computers. 
Agreement  of  the  calculator  programs  with  the  FACT  model  is 
fairly  good  when  conditions  are  within  the  design  limitations 
of  the  programs  and  environmental  conditions  are  not  unusual. 
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I.   THE  NEED  FOR  GOOD  CZ  PREDICTIONS 

A.   THE  PROBLEM  OF  THE  AIRBORNE  ASW  UNIT 

In  most  of  the  Pacific  Ocean  and  much  of  the  Atlantic 
Ocean  and  Mediterranean  Sea,  convergence  zone  (CZ)  conditions 
exist  a  majority  of  the  time,  and  they  provide  passive  acous- 
tic sensors  with  an  important  means  of  detecting  sounds  emit- 
ted from  submarine  targets .   In  some  areas  the  CZ  regions 
are  the  most  important  contact  regions  considered  in  planning 
acoustic  searches.   Obtaining  accurate  predictions  of  CZ 
sound  propagation  is  therefore  vital  to  the  success  of  acous- 
tic sensor  tactical  planners. 

Currently,  the  primary  source  of  acoustic  predictions  for 
U.  S.  Navy  units  is  the  Fleet  Numerical  Weather  Center,  Mon- 
terey, California.   Propagation  loss  profiles  for  four  stan- 
dard frequencies  and  three  source  and  receiver  depth 
combinations  are  normally  provided  in  the  ASW  Range  Predic- 
tion System  (ASRAPS)  to  air  ASW  units  when  requested.   The 
profiles,  showing  transmission  loss  (TL)  versus  distance, 
are  generated  on  a  large  digital  computer  which  uses  the 
Fast  Asymptotic  Coherent  Transmission  (FACT)  model.   This 
model  uses  as  inputs  a  linearly  segmented  sound  velocity 
profile  (SVP) ,  source  and  receiver  depths,  and  frequencies 
of  interest.   The  SVP  used  can  be  specified  by  the  user  or 
can  come  from  information  stored  at  FNWC  in  the  form  of 
historical  data.   This  stored  data  is  updated  by  bathythermo- 
graph CBT)  reports  through  a  complex  weighting  scheme  as  the 
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reports  are  received.   Without  going  into  further  detail,  it 
can  be  stated  that  predictions  produced  are  only  as  good  as 
the  data  and  the  computer  model  used,  and  only  as  timely  as 
communications  allow. 

When  entering  a  search  area,  a  problem  often  arises  con- 
cerning the  TL  profiles  obtained  from  FNWC .   Upon  taking  a 
BT  measurement,  the  unit  often  finds  the  BT  profile  used  to 
generate  the  acoustic  predictions  does  not  agree  with  the 
actual  BT  conditions  in  the  area.   If  this  situation  occurs, 
the  unit  tends  to  lose  confidence  in  the  accuracy  of  the 
predictions  and  tactical  effectiveness  is  felt  to  be  dimin- 
ished by  lack  of  good  information.   The  objective  of  this 
study  was,  therefore,  to  investigate  what  could  be  done  with 
state-of-the-art  programmable  calculators  to  improve  on  the 
in  situ  convergence  zone  predictions  available  to  air  ASW 
units . 

The  reader  only  interested  in  the  calculator  programs 
developed  may  skip  immediately  to  the  appendix. 

B.   ARE  LARGE  COMPUTERS  NECESSARY? 

In  section  5.6  of  Ref.  1,  Principles  of  Underwater  Sound, 
by  R.J.  Urick,  the  author  discusses  the  relative  merits  of 
two  theoretical  approaches  to  obtaining  wave  equation  solu- 
tions in  order  to  describe  the  distribution  of  sound  energy 
in  space  and  time.  Several  references  are  made  to  the  need 
for  digital  computers  to  produce  sound  propagation  descrip- 
tions with  either  theory.   Since  those  comments  were  made, 
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however,  there  has  been  a  revolution  in  the  capabilities  of 
small  programmable  calculators.   Although  it  is  probably  true 
that  computers  are  required  to  produce  a  complete  description 
of  sound  propagation  in  the  ocean  with  one  program,  calcula- 
tors are  capable  of  solving  the  different  modes  of  propaga- 
tion one  at  a  time  with  separate  programs  to  obtain  a  composite 
description.   Examples  of  simple  but  fairly  adequate  calculator 
programs  for  surface  duct,  bottom  bounce,  reliable  acoustic 
path,  and  deep  sound  channel  propagation  modes  are  contained 
in  Refs.  2  and  3.   These  references  also  contain  simple  models 
for  CZ  propagation,  but  they  are  based  on  a  mix  of  ray  theory, 
rule-of -thumb ,  and  empirical  data.   It  was  felt  that  a  better 
CZ  model  needed  to  be  developed. 

C.   DESIRABLE  CHARACTERISITCS  OF  A  CALCULATOR  PROGRAM  FOR 
CZ  PREDICTIONS 

1.  The  program  should  require  a  minimum  of  easily  avail- 
able input  data.   The  only  information  not  currently  available 
but  which  would  be  needed  by  an  airborne  ASW  unit  is  an  SVP 
from  the  permanent  thermocline  to  the  ocean  bottom  in  the 
unit's  search  area.   A  chart  of  this  data  could  easily  be  in- 
cluded in  the  environmental  package  carried  aboard  the  aircraft 

2.  The  program  should  be  easy  to  operate  and  not  require 
the  operator  to  have  a  great  deal  of  insight  into  the  mathe- 
matical model  or  the  internal  operations  of  the  calculator. 

3 .  The  output  should  provide  ranges  to  the  inner  and 
outer  edges  of  all  CZ  annuli  of  interest.   It  should  also  pre- 
sent the  expected  TL  for  all  frequencies  of  interest  in  each 

annulus . 
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4.  The  run  time  of  the  program  should  be  relatively 
short.   This  characteristic  recognizes  that  time  is  impor- 
tant to  an  on-station  unit. 

5.  The  program  should  be  based  on  generally  accepted 
acoustic  theory.   This  characteristic  is  desirable  because  a 
user  would  probably  have  more  confidence  in  such  a  model 
than  one  based  on  empirical  data  and  thus  applicable  only  to 

a  specific  ocean  basin.   With  empirical  models,  the  user  often 
wonders  if  the  area  he  intends  to  search  corresponds  to  the 
mean  set  of  conditions  used  to  generate  the  model  or  is  some- 
how different. 

6.  Ideally,  the  program's  performance  should  agree 
closely  with  the  generally  accepted  large  computer  models 
currently  in  use. 
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II.   USING  ICAPS  AS  THE  STANDARD  FOR  COMPARISON 

A.  RATIONALE  FOR  ONLY  ONE  LOCATION  PER  OCEAN  BASIN 

As  will  be  demonstrated,  the  CZ  characteristics  of  the 
three  locations  studied  vary  considerably.   The  deep  sound 
channels  which  produce  CZ  phenomena  are  quite  different  as 
are  the  ranges  from  source  to  CZ  annuli.   The  objective  of 
this  study  was  to  produce  a  mathematical  model  for  CZ  pre- 
dictions for  use  on  small  programmable  calculators.   It  was 
reasoned  that  if  the  model  would  work  for  the  different  con- 
ditions of  the  three  locations  studied,  it  would  work  for  all 
of  the  variations  to  be  expected  within  any  one  of  the  ocean 
areas. 

B.  DESCRIPTION  OF  ICAPS 

The  Integrated  Carrier  ASW  Prediction  System  (ICAPS)  is 
a  passive  and  active  acoustic  prediction  system  developed 
for  installation  aboard  aircraft  carriers  and  other  large 
naval  vessels  which  have  digital  computers.   It  contains 
four  sets  of  historical  environmental  data,  one  each  for  the 
North  Pacific,  North  Atlantic,  and  Indian  Oceans,  and  one 
for  the  Mediterranean  and  Black  Seas.   It  also  contains  sev- 
eral production  programs  for  predicting  naval  sonar  system 
performance.   The  FACT  model  is  used  in  the  passive  sensor 
predictions.   This  is  the  same  model  used  at  FNWC  for  ASRAPS. 
Reference  4  contains  a  description  of  the  installation  and 
operation  of  ICAPS  in  the  IBM  3  60  Computer  Center  at  the 
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Naval  Postgraduate  School,   Reference  5  contains  a  descrip- 
tion of  the  mathematics  used  in  the  FACT  model. 

C.   DESCRIPTION  OF  HISTORICAL  EN VI RONMENT AL  DATA  FILES 

Figures  1(a)  through  1CD  depict  twelve  sound  velocity 
profiles  produced  by  ICAPS  from  its  historical  environmental 
data  files.   Figures  ICa)  through  1(d)  show  SVP  information 
for  the  months  of  February,  May,  August,  and  November  for 
40N  140W  in  the  Pacific  Ocean.   Figures  1(e)  through  1(h) 
are  the  same  information  for  3 IN  69W  in  the  Atlantic  Ocean, 
and  likewise,  Figures  l(i)  through  111)  are  for  36N  18E  in 
the  Mediterranean  Sea. 

The  historical  data  files  used  to  produce  these  profiles 
consist  of  temperature  and  salinity  values  for  over  thirty 
depths,  for  four  seasons  of  the  year,  and  for  many  locations 
spaced  at  one  to  five  degree  latitude  and  longitude  inter- 
vals in  each  ocean  area  covered.   When  specific  latitude, 
longitude,  and  date  are  specified,  interpolations  are  per- 
formed to  produce  the  approximate  temperature  and  salinity 
profiles  to  be  expected  at  that  location  and  date.   This  in- 
formation is  then  converted  to  an  SVP.   The  output  from  this 
portion  of  the  system  consists  of  seven  columns  of  values, 
one  each  for  depth  in  meters  and  feet,  temperature  in  Cel- 
sius and  Fahrenheit,  salinity,  and  sound  velocity  in  meters 
per  second  and  feet  per  second.   The  depths  associated  with 
these  quantities  begin  at  ten  meter  intervals  near  the  ocean 
surface  and  gradually  increase  through  25,  50,  100,  250,  500, 
and  1,000  meter  intervals  as  depth  increases.   The  last  line 
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of  values  is  for  the  ocean  bottom  depth  which  was  part  of 
the  input  data. 

D.   COMPARISON  OF  DEEP  SOUND  CHANNEL  CHARACTERISTICS 

Figure  2  depicts  the  deep  sound  channel  (DSC)  portion 
of  the  May  SVP  for  the  Pacific,  Atlantic  and  Mediterranean 
coordinates  mentioned  earlier  in  a  composite  graph  drawn  to 
scale.   As  can  be  seen  in  that  figure,  DSC  characteristics 
of  the  three  areas  differ  considerably.   The  vertical  ex- 
tent of  the  channels  varies  from  1100  meters  in  the  Mediter- 
ranean to  4200  meters  in  the  Atlantic.   The  change  in  sound 
velocity  between  sonic  layer  depth  (SLD)  and  DSC  axis  (point 
of  minimum  velocity)  varies  from  15  m/sec  in  the  Mediterra- 
nean to  38  m/sec  in  the  Atlantic.   The  depth  of  the  DSC  axis 
varies  from  100  meters  in  the  Mediterranean  to  1300  meters 
in  the  Atlantic.   Pacific  Ocean  values  are  between  the  others 
for  all  of  those  characteristics.   Sound  velocity  near  the 
surface  is  much  greater  in  the  Atlantic  and  Mediterranean 
than  in  the  Pacific,  and  sound  velocity  near  the  bottom  of 
the  three  basins  (not  shown  in  the  figure)  is  about  6  m/sec 
greater  in  the  Atlantic  than  in  the  Pacific  and  about  52 
m/sec  greater  in  the  Mediterranean  than  in  the  Pacific  at 
equal  depths.   Also  note  the  subsurface  sound  channel  lo- 
cated about  100  to  500  meters  below  the  surface  in  the  At- 
lantic profile. 
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E.   COMPARISON  OF  CZ  CHARACTERISTICS  IN  THREE  OCEANS 
1 .   Method  Used  in  Obtaining  Data  for  Comparison 
There  were  two  primary  objectives  in  gathering 
twelve  ICAPS  runs  from  each  of  the  three  ocean  areas.  First, 
it  was  desired  to  obtain  sufficient  data  to  determine  which 
CZ  characteristics  are  common  to  all  areas  and  which  charac- 
teristics are  peculiar  to  specific  basins.   Secondly,  it  was 
desired  to  obtain  a  standard  of  comparison  for  any  calculator 
program  which  might  be  developed.   To  fulfill  the  first  ob- 
jective, it  was  decided  to  keep  the  input  variables  the  same 
in  all  areas,  varying  them  one  at  a  time,  in  order  to  better 
compare  the  differences  observed  in  the  various  runs.   For 
each  of  the  three  locations ,  the  inputs  varied  were  season 
of  the  year  and  source  and  receiver  depth  combination.   Re- 
ceiver depths  of  60  and  300  feet  and  source  depths  of  60  and 
4Q0  feet  were  used.   Each  of  the  ICAPS  outputs  consisted  of 
TL  profiles  for  four  frequencies  out  to  a  range  of  250  kyds . 

Originally,  it  was  intended  to  collect  twelve  data 
from  each  profile.   These  data  were  to  be  the  range,  width, 
CZ  gain,  and  Transmission  Loss  for  each  of  the  first  three 
convergence  zone  annuli .   As  it  turned  out,  somewhat  less 
data  was  collected  and  tabulated.   There  were  several  rea- 
sons for  this.   First,  the  February  SVP  in  the  Mediterranean 
contained  no  sound  channel  and  therefore  no  convergence  zones 
existed.   Secondly,  all  of  the  third  CZ  data  for  the  Atlan- 
tic was  thrown  out  on  the  grounds  that  it  was  almost  always 
the  same  and  that  it  was  inconsistent  with  information  from 
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the  first  two  CZ  annuli  in  any  particular  profile.   The 
reason  for  this  occurrence  is  not  know.   Finally,  it  was 
impossible  to  obtain  some  of  the  desired  data  because  of 
the  smooth  way  in  which  the  CZ  path  blended  with  other  com- 
petitive propagation  modes.   One  could  not  tell  what  was 
CZ  and  what  was  not  in  those  cases. 
2 .   CZ  Range  and  Width  Analysis 

The  transmission  loss  profiles  produced  by  ICAPS 
are  presented  in  two  formats,  a  table  of  TL  values  for  each 
kiloyard  of  range  from  the  source  and  a  graph  of  the  same 
information.   Because  the  TL  values  are  tabulated  at  kilo- 
yard  intervals,  it  is  impossible  to  be  more  accurate  than 
that  interval  in  determining  where  a  CZ  begins  and  ends . 
Also,  it  was  difficult  to  be  consistent  in  picking  the 
points  representing  the  edges  of  CZ  annuli  because  of  the 
variety  of  graph  shapes,  TL  levels,  and  other  propagation 
mode  interferences.   In  any  event,  an  attempt  was  made  to 
satisfy  one  basic  criterion  in  choosing  leading  and  trailing 
edges  of  the  annuli:   Do  these  ranges  best  represent  the 
apparent  location  of  the  annulus  regardless  of  the  TL  levels 
involved?   Admittedly,  the  ranges  picked  were  often  based 
on  subjective  judgement,  and  it  cannot  be  stated  with  com- 
plete certainty  that  only  CZ  mode  propagation  contributed 
to  the  TL  peaks  judged  to  be  the  CZ  annuli. 

Table  I  contains  the  CZ  range  and  width  data  that 
could  be  gleaned  from  the  ICAPS  profiles.   In  the  table, 
RCZi  is  the  range  to  the  inner  edge  of  the  first,  second,  or 
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third  CZ  annulus  to  the  nearest  one  half  nautical  mile.   CZW 
is  the  width  of  the  annulus  also  the  nearest  one  half  nauti- 
cal mile.   The  third  column  of  numbers  is  the  ratio  of  CZW 
to  the  range  of  the  outer  edge  of  the  annulus  (RCZo) ,  ex- 
pressed as  a  percentage. 

After  carefully  studying  this  data,  the  following 
conclusions  were  made  concerning  CZ  propagation: 

(1)  The  range  to  the  first  CZ  is  approximately 
14  to  18  nm  at  the  Mediterranean  location,  23  to  27  nm  at 

the  Pacific  location,  and  33  to  35  nm  at  the  Atlantic  location 

(2)  Range  to  the  CZ  decreases  and  annulus 
width  increases  as  source  and  receiver  get  deeper  in  all 
cases . 

(3)  The  ranges  to  the  second  and  third  annuli 
are  approximately  whole  number  multiples  of  the  ranges  to  the 
inner  and  outer  edges  of  the  first  annulus  in  all  cases. 

(4)  The  range  or  width  of  a  CZ  annulus  does 
not  appear  to  have  any  significant  frequency  dependence. 

3 .   CZ  Gain  and  Transmission  Loss  Analysis 

Convergence  zone  gain  is  defined  as  the  difference 
between  the  transmission  loss  expected  under  conditions  of 
spherical  propagation  and  the  actual  transmission  loss  ob- 
served.  This  definition  is  expressed  in  Eq.  (1). 

G  =  20  log(r)  +  a(r)  -  TL  CD 

In  this  equation,  G  is  the  CZ  gain,  r  is  the  range  to  the 
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CZ  annulus ,  a  is  the  attenuation  coefficient  associated  with 
the  frequency  of  interest,  and  TL  is  the  actual  transmission 
loss  observed  in  the  CZ  annulus  for  that  frequency.   All 
terms  in  Eq.  (1)  are  in  decibels  (dB) . 

In  actual  convergence  zones,  TL  (and  therefore  gain) 
is  by  no  means  a  constant  value.   Contributions  of  several 
possible  propagation  paths  at  any  one  point  and  the  time 
varying  nature  of  sound  paths  in  the  ocean  cause  coherence 
effects  to  exist.   These  effects  make  TL  vary  in  both  space 
and  time.   Coherence  effects  are  more  pronounced  at  lower 
frequencies  (longer  wavelengths)  where  the  time  varying  ef- 
fects are  small  compared  to  spatially  distributed  effects. 
In  ICAPS ,  the  more  predictable  coherence  conditions  are  in- 
cluded in  the  mathematical  model. 

Since  a  single  TL  value  was  desired  for  the  envi- 
sioned calculator  model,  an  attempt  was  made  to  pick  the 
"average"  TL  in  the  ICAPS  CZ  annuli.   As  with  the  range  es- 
timates, this  called  for  subjective  judgement.   Figure  3 
shows  a  typical  ICAPS  CZ  presentation  which  has  coherence 
effects  in  evidence.   The  figure  suggests  how  an  "average" 
TL  was  chosen  as  best  representing  that  annulus.   Two  levels 
were  chosen  (labeled  high  and  low  in  the  figure)  which  bracket 
the  majority  of  the  TL  points  within  the  annulus.   The  ap- 
proximate midpoint  between  those  levels  was  then  picked  as 
"the"  TL  for  that  CZ . 

As  an  extra  point  of  interest,  the  high  and  low  TL 
levels  were  studied.   It  was  noted  that  ICAPS  predicts  TL 
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Ran^e  (run) 


Figure  3-   Estimating  Transmission  Loss  in  a  CZ 
Annulus  from  an  ICAP3  TL  Frofile. 
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variations  from  about  i  lOdB  around  the  "average"  level.   If 
this  is  truly  representative  of  CZ  coherence  effects,  an  ASW 
unit  armed  only  with  an  estimate  of  the  "average"  TL  in  a 
certain  CZ  annulus  should  expect  to  see  variations  of  about 
that  magnitude  around  the  estimate  in  hand. 

Using  TL  levels  estimated  by  the  procedure  described 
above,  and  employing  Eq.  CD /  CZ  gain  values  predicted  by 
ICAPS  were  obtained  and  tabulated.   The  values  produced  are 
contained  in  Tables  II  and  III.   Table  II  shows  all  of  the 
data  from  the  Pacific  location.   Table  III  contains  only  300 
Hz  data  from  the  Atlantic  and  Mediterranean  locations.   (50, 
850,  and  1700  Hz  data  were  omitted  from  Table  III  because  it 
became  obvious  during  data  collection  that  G  is  not  frequency 
dependent. ) 

Again  after  careful  study,  the  following  conclusions 
were  drawn  concerning  CZ  gain: 

(1)  CZ  gain  values  range  from  about  eight  to 
twenty  dB  in  all  three  areas  observed. 

(2)  CZ  gain  is  the  same  value  for  first, 
second,  and  third  CZ  in  any  given  case. 

(3)  In  general,  CZ  gain  is  independent  of 
frequency.   An  exception  to  this  conclusion  is  that  at  low 
frequency  (below  300Hz) ,  especially  when  source  or  receiver 
or  both  are  above  the  SLD  and/or  near  the  surface,  there  is 
apparently  somewhat  less  gain  than  evident  for  higher  fre- 
quencies.  This  difference  is  probably  due  to  stronger  dif- 
fraction of  the  longer  wavelengths. 
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C4)   CZ  gain  seems  to  be  highest  when  source 
and  receiver  are  at  or  near  the  same  depth. 
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Month     Rcvr/Tgt 

Freq 
(Hz) 

CZ  Gain  (dB) 
1st  CZ     2nd  CZ 

SLD(ft)      (ft) 

3rd  CZ 

50 

7 

7 

6 

60/60 

300 
850 

14 
14 

13 
13 

12 
13 

1700 

15 

15 

12 

50 

9 

9 

9 

246        6°/400 

300 
850 

12 

11 

12 

11 

11 
12 

1700 

11 

11 

12 

50 

12 

11 

10 

300/400 

300 

850 

12 

10 

12 
11 

11 
12 

1700 

12 

10 

11 

50 

16 

17 

16 

60/60 

300 
850 

16 
16 

17 
17 

16 
16 

1700 

16 

17 

17 

50 

10 

10 

12 

^y|        60/400 

300 
850 

10 
11 

11 
12 

13 
15 

1700 

9 

11 

12 

50 

14 

14 

14 

300/400 

300 
850 

13 
13 

13 

14 

14 
11 

1700 

12 

11 

13 

Table  II.   ICAPS  CZ  Gain 

Data 

Observed  at  the  Pacif 

ic  Ocean 

Location.       (Page    1   of    2 
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Month     Rcvr/Tgt 
SLD(ft)     (ft) 


Freq 

CZ 

Gain  (dB) 

(Hz) 

1st  CZ 

2nd  CZ 

3rd  CZ 

50 

11 

11 

13 

300 

14 

13 

12 

850 

15 

18 

20 

1700 

16 

17 

18 

60/60 


50  7  7  8 

AUG                   300  9  12  13 

—        60/40°      850  12  14  15 

1700  12  13  15 


50  14  14  14 

300/400     300  11  14  12 

850  11  12  15 

1700  12  13  11 


60/60 


50  8  7  6 

300  9  12  12 

850  15  18  16 

1700  14  15  15 


50  9  9  8 

NOV        -n  ,,nn      300  12  13  11 

"T8        60/40°      850  11  11  11 

1700  11  11  11 


50  11  11  11 

300/400     30°  12  X1  1X 

JUU/4UU      85Q  9  ]_1  12 

1700  10  9  11 


Table  II.   ICAPS  CZ  Gain  Data  Observed  at  the  Pacific  Ocean 
Location.   (Page  2  of  2) 
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Ocean    Month 

Rcvr/Tgt 

300 

Hz  CZ 

Gain 

SLD(ft) 

(ft) 

1st  CZ    2 

nd  CZ 

3rd  CZ 

FEB 

60/60 

17 

14 

328 

60/400 

13 

13 

300/400 

17 

18 

MAY 
0 

60/60 

19 

14 

60/400 

16 

12 

u 

300/400 

13 

17 

E-» 

J 

AUG 
0 

60/60 

15 

17 

% 

60/400 

10 

13 

300/400 

14 

15 

NOV 
0 

60/90 

13 

11 

60/400 

12 

10 

300/400 

13 

12 

MAY 

60/60 

14 

13 

16 

10 

60/400 

12 

10 

10 

300/400 

17 

18 

20 

z 

< 

AUG 

60/60 

17 

19 

17 

z 

Pi 

o 

60/400 

11 

9 

10 

300/400 

w 

Eh 

H 

Q 

NOV 

60/60 

11 

13 

15 

S 

10 

60/400 

10 

11 

12 

300/400 

12 

14 

15 

Table  i: 

CI.   ICAPS  300Hz 

CZ  Gain  Data  Observed 

at  Atlantic 

and  Mediterranean  Locations 
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III.   CZ  RAY  THEORY  ANALYSIS  AND  MODEL  DEVELOPMENT 

A.   CZ  RANGE  AND  WIDTH 

It  was  decided  to  use  ray  tracing  as  the  method  for  de- 
termining CZ  range  and  width  because  of  the  simplicity  of 
the  mathematics  involved  and  because  of  the  intuitive  appeal 
of  sound  rays  depicting  the  propagation  of  sound.   The  alter- 
native approach,  that  of  normal  mode  theory,  was  rejected  on 
the  grounds  that  it  would  be  much  more  complicated,  requiring 
capabilities  far  beyond  those  available  in  the  calculators 
at  hand. 

Figure  4  shows  four  sound  rays  of  particular  interest  in 
CZ  propagation.   The  order  of  these  rays  is  described  for  the 
"typical"  case  in  the  following  discussion.   An  "atypical" 
case  will  be  mentioned  later. 

Ray  #1  departs  the  SLD  at  zero  degree  depression  angle. 
It  reaches  its  greatest  depth  at  the  bottom  of  the  DSC  and 
returns  to  the  SLD  at  some  particular  range  and  at  horizontal 
incidence.   The  horizontal  range  from  SLD  to  SLD  is  termed 
cycle  distance.   The  cycle  distance  for  this  ray  is  desig- 
nated r0  . 

Ray  #2  is  the  next  ray  of  interest  found  as  the  departure 
angle  from  the  SLD  is  increased  downward.   This  ray  passes 
down  through  and  below  the  bottom  of  the  DSC  before  turning 
back  upward.   It  returns  to  the  SLD  at  the  shortest  range 
from  the  starting  point  of  any  ray  within  the  bundle  of  rays 
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undergoing  CZ  refraction.   This  range  is  designated  r  .  . 

3  •  min 

The  angle  of  departure  for  this  ray  is  designated  0   .  . 
Each  ray  between  #1  and  #2  crosses  all  of  the  previous  (les- 
ser departure  angle)  rays  on  its  way  up  from  its  lowest 
depth. 

As  departure  angle  from  the  SLD  is  further  increased, 
the  next  sound  ray  of  interest,  #3,  is  located.   This  ray 
has  a  cycle  distance  equal  to  r„ .   Its  maximum  depth  is 
greater  than  that  for  ray  #2.   It  departs  from  and  arrives 

back  at  the  SLD  at  an  angle  designated  9     .   Rays  between 

rswp      J 

#2  and  #3  do  not  cross  each  other,  but  they  do  cross  the 
earlier  rays  on  their  way  back  up  to  the  SLD. 

In  the  CZ  annulus ,  the  rays  between  #1  and  #2  sweep  in- 
ward toward  the  source  as  departure  angle  increases.   After 
ray  #2  they  sweep  out  away  from  the  source  as  angle  increases 
further.   For  this  reason,  the  region  formed  by  rays  between 
#1  and  #3  is  called  the  reswept  zone. 

Finally,  as  angle  of  departure  from  the  SLD  is  increased 
to  maximum  angle  for  CZ  propagation,  we  observe  ray  #4.  This 
ray  turns  upward  at  a  depth  equal  to  the  water  column  depth 
at  that  location.   It  returns  to  the  SLD  at  the  greatest  dis- 
tance of  all  CZ  refracted  rays.   Rays  departing  the  SLD  at 
angles  greater  than  that  for  ray  #4  would  be  reflected  off 
the  bottom  and  are  not  of  significance  for  the  CZ  propagation 
path. 

As  mentioned  earlier,  this  progression  of  rays  exists 
in  a  "typical"  CZ  situation.   If,  however,  the  ocean  bottom 
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were  more  shallow,  cycle  distance  for  ray  #4  would  be  re- 
duced.  If  the  bottom  were  shallow  enough,  ray  #4's  cycle 
distance  would  be  less  than  r  .   In  that  case,  the  reswept 
zone  would  be  reduced  to  the  region  between  rays  #2  and  #4. 
This  situation  is  called  the  "atypical"  case. 

B.   RANGE  AND  WIDTH  MODEL 

In  the  calculator  programs  developed,  provision  is  made 
for  entering  and  storing  a  five  point  sound  velocity  profile 
which  defines  the  DSC  only.   The  first  depth  and  velocity 
pair  entered  (D,,C,)  equate  to  the  appropriate  values  found 
at  the  SLD.   The  fifth  depth  entered  (D_)  is  the  depth  at 
the  bottom  of  the  DSC  where  sound  velocity  is  equal  to  that 
at  the  SLD.   The  other  three  depth/velocity  pairs  must  be 
picked  subjectively  from  a  graph  of  the  SVP  of  interest.   If 
a  mixed  layer  exists,  the  gradient  in  that  layer  is  taken  to 
be  0.Q2  sec    (a  purely  pressure  induced  gradient  to  two 
place  accuracy) .   The  program  calculates  the  four  layer  gra- 
dients within  the  DSC  profile  entered,  and  uses  the  fourth 
(deepest)  layer  gradient  in  ray  calculations  that  occur  below 
the  DSC.  It  would  have  been  desirable  to  allow  several  more 
points  in  the  SVP,  but  calculator  data  storage  capacity  and 
program  step  limitations  preclude  more  than  five  depth/veloc- 
ity pairs. 

The  overall  scheme  used  to  predict  CZ  range  and  width  is 
to  trace  a  series  of  rays  starting  at  the  SLD  with  a  zero  de- 
pression angle  ray.   That  first  ray  yields  r  which  is  stored 
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Then  an  iterative  process  is  begun  in  which  the  angle  is  in- 
cremented and  each  succeeding  cycle  distance  determined  is 

compared  to  the  previous  one  until  r  .   and  9   .   are  found. 

mm       rnun 

(0   .   is  stored  for  use  in  the  CZ  gain  and  TL  program,  to 
rmm  3  c      r    ' 

be  discussed  later.)   Corrections  are  then  made  to  r  .   and 

mm 

rQ  to  account  for  surface  duct  effects  (if  any)  and  source 

and  receiver  depth  separation  from  the  SLD.   Range  to  the 

inner  edge  of  the  CZ  is  r  .   plus  corrections,  and  range  to 
r  min  r  ^ 

the  outer  edge  is  r~  plus  corrections. 

In  the  first  attempt  to  produce  a  calculator  program, 
the  cycle  distance  of  ray  #4  (the  ray  just  grazing  the  bot- 
tom) was  compared  to  r_ .   The  greater  of  the  two  was  picked 
as  the  basic  distance  for  determining  range  to  the  outer 
edge  of  a  CZ.   Later  on,  this  portion  of  the  program  had  to 
be  deleted  to  save  program  steps.   The  final  programs  de- 
veloped ignore  bottom  depth  and  do  not  include  rays  outside 
the  reswept  zone  in  determining  annular  width.   This  is  prob- 
ably a  shortcoming  of  the  programs  but  the  seriousness  of 
the  errors  it  causes  will  not  be  known  without  further  study. 

A  commonly  applied  rule-of-thumb  states  there  must  be 
a  minimum  30  0  fathoms  of  depth  excess  (water  column  below 
the  DSC)  in  order  to  have  "reliable"  CZ  conditions.   It  was 
observed  that  a  fully  developed  reswept  zone  existed  in 
every  case  in  the  locations  studied,  and  separate  calcula- 
tions showed  that  somewhat  less  than  30  0  fathoms  depth  ex- 
cess was  required  to  complete  the  zone.   Therefore,  a  program 
user  should  consider  the  300  fathom  rule-of-thumb  before 
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running  the  range  and  width  program.   With  less  than  300 
fathoms  depth  excess,  the  possibility  exists  for  an  "atypical" 
CZ  propagation  situation  where  the  reswept  zone  is  reduced  in 
width  due  to  bottom  ray  limiting. 

Another  program  shortcoming  involves  an  assumption  that 
both  source  and  receiver  would  be  more  shallow  than  the  sec- 
ond DSC  SVP  point  chosen  (depth  D~)  .   In  other  words,  the 
programs  were  designed  to  allow  for  source/receiver  depths 
within  the  mixed  layer  or  the  first  isogradient  layer  below 
the  SLD.   After  the  five  point  SVP  is  entered  and  the  grad- 
ients computed,  source  and  receiver  depths  are  entered  and 
converted  to  velocities.   The  programs  determine  these  veloc- 
ities (C   and  C  )  by  subtracting  an  appropriate  amount  from 
the  velocity  at  the  SLD.   The  amount  subtracted  is  determined 
by  depth  separation  from  the  SLD  and  by  the  gradient  in  either 
the  ML  or  the  first  layer  below  the  SLD.   If  source  or  re- 
ceiver depth  is  greater  than  D_ ,  sound  velocity  should  be  de- 
termined by  correcting  C-  (the  velocity  at  D-)  and  by  using 
g~  (the  second  layer  gradient) .   Since  this  is  not  done,  ve- 
locities for  source/receiver  depths  below  D~  will  be  in  error 
(usually  too  low) .   Source  and  receiver  velocity  errors  are 
carried  over  into  Arc  and  Ar   range  corrections.   If  the 
velocities  are  too  low,  the  range  corrections  will  be  too 
large.   This  is  normally  a  rather  insignificant  source  of 
total  range  error,  however,  since  Ars  and  Ar   errors  will  be 
a  small  fraction  of  the  magnitude  of  those  terms  and  because 
the  range  correction  terms  are  small  to  begin  with. 


The  mathematics  of  ray  tracing  in  isogradient  layers  is 
quite  straightforward.   By  Snell's  Law, 


'2    _   (A  constant 


(2) 


cos  9.    cos  0_    for  each  ray) 

the  angle  of  a  ray  departing  a  layer  can  be  determined  from 
the  angle  of  entry  into  that  layer.   In  Eq.  (2)  ,  C-,  is  the 
sound  velocity  where  the  sound  ray  enters  a  layer,  9,  is  the 
angle  of  entry,  C_  is  the  sound  velocity  where  the  ray  de- 
parts the  layer,  and  9_  is  the  angle  of  departure. 

Rays  travel  in  circular  arcs  within  constant  gradient 
layers,  and  the  radius  of  curvature  is: 


R  = 


g.  cos 


1 


(3) 


C,  and  9,  are  as  defined  above  and  g,  is  the  gradient  within 
the  layer  (in  this  case,  layer  1)  . 

Finally,  the  horizontal  distance  traveled  by  a  ray 
while  traversing  a  layer  is: 


Ar  = 


R  (sin  9  2  -  sin  9 , ] 


(4 


Figure  5  demonstrates  an  example  application  of  Eqs . 
2  through  4.   It  should  be  noted  that  absolute  value  signs 
are  used  in  Eqs.  3  and  4  because  the  gradient  in  Eq.  3  and 
the  difference  of  sines  in  Eq.  4  may  be  positive  or  negative, 
while  R  and  Ar  are  always  positive. 

The  programs  use  these  equations  to  compute  the  horizon- 
tal range  increments  each  ray  accumulates  within  the  four 
layers,  doubles  each  term  (.to  account  for  the  downward  and 
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upward  passes  through  each  layer) ,  and  then  sums  the  terms 
to  obtain  cycle  distances.   The  fourth  layer  requires  a 
slightly  different  treatment  because  the  rays  become  horizon- 
tal and  then  turn  back  upward  within  that  layer.   Essentially 
the  same  formulas  are  used,  however.   The  equations  are  also 
used  to  compute  the  range  correction  terms . 

In  considering  the  various  possible  ray  paths  between 
source  and  receiver,  it  was  decided  there  were  four  basic 
situations  which  could  occur: 

1)  No  mixed  layer,  both  source  and  receiver 
below  the  SLD .   (Deep/Deep) 

2)  Mixed  layer  present,  both  source  and  receiver 
below  the  SLD.   (Deep/Deep/ML) 

3)  Mixed  layer  present,  both  source  and  receiver 
within  the  layer.   (Shal/Shal) 

4)  Mixed  layer  present,  source  or  receiver  above 
the  SLD,  the  other  below.   (Crosslayer) 

The  only  difference  between  the  first  two  cases  is  the  mixed 
layer  effect  in  case  2.   That  effect  causes  a  widening  of 
annuli  due  to  spreading  of  sound  rays  as  they  travel  up  to 
the  surface  and  back  down  to  the  SLD  within  the  layer.   The 
mixed  layer  effect  is  also  included  in  the  third  and  fourth 
cases  above. 

It  was  originally  intended  to  include  all  four  cases  in 
one  range  prediction  program.   Again  due  to  calculator  limi- 
tations, it  was  necessary  to  use  two  programs  to  cover  the 
four  possibilities.   The  first  range  program  (labeled  Deep/ 
Deep)  is  for  cases  1)  and  2)  above  when  both  source  and  re- 
ceiver are  below  the  SLD  whether  or  not  an  ML  exists.   The 


55 


second  range  program  (labeled  Shal/Shal  or  Crosslayer)  is 
for  use  in  cases  3)  and  4)  above  when  source  or  receiver  or 
both  are  above  the  SLD. 

Formulas  used  to  determine  range  to  inner  edge  of  the 
first  CZ  CRCZi)  and  range  to  the  outer  edge  of  the  first  CZ 
(RCZo)  follow: 

RCZi  =  rmin  "  ArS  "  ArR  Deep/Deep 

=  r  .   +  Arc  +  Ar_  Shal/Shal 

mm      S      R  ' 

=  r  .   +  Ar   -  Ar_  Crosslayer 

mm     s     R  2 

RCZo  =  rQ  +  (2  Ar  }  +  Arg  +  ArR       Deep/Deep 

0  Deep/Deep/ML 

=  rQ  +  2  Ar   +  Ar   -  Ar  Shal/Shal 

=  r_  +  2  Ar.  +  Ar   +  Ar  Crosslayer 

In  these  equations,  r     and  r,  have  been  previously  defined, 
n  mm       0  e  jt  i 

Ar  and  Ar   are  the  respective  horizontal  range  corrections 
which  account  for  source  and  receiver  depth  separation  from 
the  SLD,  and  2  Arn  is  the  correction  for  mixed  layer  effect. 
Figures  6(a)  through  6(f)  (not  to  scale)  depict  the  RCZi  and 
RCZo  formulas  in  graphic  form.   Ranges  to  second  and  subse- 
quent CZ  annuli  are  taken  to  be  integer  multiples  of  the 
ranges  produced. 

Two  items  of  interest,  both  evident  in  Figs.  6(a)  -  6(f) , 
are  worth  mentioning  at  this  point.   First,  acoustical  reci- 
procity is  envoked  and  the  more  shallow  of  source  and  re- 
ceiver is  always  treated  as  "source"  of  the  sound  rays  within 
the  calculator  programs.   Secondly,  only  those  sound  rays 
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which  experience  no  more  than  one  ocean  surface  reflection 
between  source  and  receiver  are  considered  in  this  model. 
Both  of  these  conventions  are  commonly  applied  to  ray  trac- 
ing models.   Although  they  theoretically  have  little  or  no 
effect  on  model  results,  they  greatly  simplify  the  work  of 
programming  a  ray  tracing  model. 

C.   CZ  GAIN  AND  TRANSMISSION  LOSS  MODEL 

In  general,  transmission  loss  is  defined  as  ten  times 
the  logarithm  of  the  ratio  of  sound  intensities  measured  at 
one  meter  from  a  source  and  at  range  r  from  that  source. 

TL  =  10  log  —=■ 
r 

Intensity  has  units  of  power  per  unit  area.   The  change  in 
intensity  between  one  meter  and  range  r  is  due  to  geometric 
spreading  of  the  power  over  a  different  amount  of  area  and 
due  to  attenuation  of  some  of  the  power  through  absorption, 
scattering,  diffusion,  etc. 

In  ray  tracing  theory,  it  is  assumed  there  is  no  sound 
power  transfer  across  sound  rays.   Therefore,  the  power  flow- 
ing from  a  source  between  two  sound  rays  remains  between 
those  rays  and  travels  out  in  a  direction  parallel  to  the 
ray  paths.   Determining  the  portion  of  transmission  loss  due 
to  geometric  spreading  (TL  )  under  this  assumption  reduces  to 
finding  ten  times  the  logarithm  of  the  ratio  of  areas  (at 
range  r  and  at  one  meter)  penetrated  by  the  power  between 
the  two  rays  perpendicular  to  the  direction  of  travel. 
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■  TLg  "  10  io?  5^ 

A  mathematical  development  of  this  technique  is  contained  on 
pages  119-121  of  Ref.  1. 

Figure  7  shows  how  this  method  was  adapted  for  use  in 
the  CZ  Gain  and  Transmission  Loss  portion  of  the  calculator 
model  developed.   The  area  CA, )  at  one  meter  from  the  source 
is  the  product  of  area  height  and  area  circumference.   The 
sound  rays  bounding  the  area  above  and  below  are  the  minimum 
and  maximum  departure  angle  rays  which  produce  the  reswept 
zone  in  the  CZ  annulus .   The  angular  spread  of  those  rays 
C.A8)  in  radian  units  times  the  sphere  radius  (1  meter)  is 
the  area  height.   Cosine  of  the  average  angle  of  departure 
of  the  rays  (9,)  times  the  sphere  radius  times  2tt  is  circum- 
ference of  the  area.   Therefore: 

A,  =  2tt  Ad  cos  9   (m  ) 

In  the  CZ,  the  area  (A~)  over  which  the  same  power  is 
distributed  is  also  found  by  a  product  of  area  circumference 
and  area  width.   Circumference  is  2 it  times  range  to  the  CZ 
(RCZi) .   Width  of  the  area  perpendicular  to  the  sound  rays 
is  the  product  of  CZ  annulus  width  (CZW)  and  the  sine  of  the 
average  angle  of  arrival  of  the  rays  at  the  receiver  depth 
(9  ) .   Therefore: 

2 
A?  =  2tt  RCZi  CZW  sin  9  2  Cm  ) 

and  the  geometric  TL  expression  becomes: 
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RCZi  CZW  sin  6 

TL   =  10  log  7-3 5 - 

g        3     A9  cos  9, 

Substituting  this  expression  back  into  Eq.  1,  which  is 

the  definition  of  CZ  gain,  and  reducing  to  simplest  form 

yields  the  algorithm  used  to  determine  G  in  the  calculator 

model : 

RCZi  A6  cos  0, 

G  =  10   log  r^rrz . 5 (b) 

*    CZW  sin  9 

To  implement  this  algorithm,  the  program  has  only  to 
determine  the  angular  terms  since  RCZi  and  CZW  are  available 
from  the  range  program  results .   After  one  of  the  range  pro- 
grams has  been  run,  the  user  loads  the  G  and  TL  program  into 
calculator  memory  without  altering  the  contents  of  the  data 
storage  registers  left  from  the  range  program.   Then  the 
iterative  ray  tracing  process  begun  in  the  range  program  is 
continued  in  the  gain  program.   The  angle  of  departure  of 

sound  rays  from  the  SLD  is  incremented  beyond  0      (left  in 
1  J  rmin 

storage  from  the  range  program)  and  cycle  distances  produced 

are  checked  for  approximate  equality  with  r» .   In  this  way, 

the  ray  which  completes  the  reswept  zone  is  found,  and  its 

angle  of  departure  from  the  SLD  is  9     .   Then  9      and 
3         r  rswp         rswp 

zero  degrees  (the  angle  for  the  ray  producing  cycle  distance 

rn)  are  converted  to  angles  of  departure  from  the  source 

depth,  9C_.  and  9C^  respectively,  and  angles  of  arrival  at 
SR        oO 

the  receiver  depth,  9T5_,  and  9_._  respectively,  using  Snell's 

RR        RO 

law.   The  angular  terms  in  Eq.  5  are  then  computed  using 
the  following  formulas: 
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A9    9SR  +  9S0         Deep  "source" 


Shal  "source" 
Deep  "source" 


SR 

SO 

9SR 

- 

6  so 

2 

6SR 

+ 

9so 

2 

1 

8R0 

+ 

9RR 

Shal  "source" 


Recall  that  "source"  in  the  model  refers  to  the  more  shallow 
of  source  and  receiver.   Therefore,  the  deep  "source"  forms 
of  these  formulas  are  used  only  after  using  the  Deep/Deep 
range  program.   In  all  other  cases,  the  "source"  is  consid- 
ered to  be  shallow.   Figure  8  depicts  these  angular  relation- 
ships for  the  various  depth  conditions. 

It  should  be  pointed  out  there  are  two  inherent  errors 
in  the  angular  quantities  determined.   First,  the  possible 
source  and  receiver  sound  velocity  errors  mentioned  earlier 
could  cause  the  gain  algorithm  angles  to  be  slightly  off. 
This  would  only  occur  if  depths  greater  than  D-  were  entered 

for  source  or  receiver  or  both.   Secondly,  9      is  found 

1         rswp 

for  the  ray  which  has  cycle  distance  equal  to  rn  at  the  SLD. 
Since  the  actual  CZ  ray  bundle  departs  from  the  source  depth 
(vice  SLD)  and  arrives  at  the  receiver  depth  (vice  SLD) ,  the 
ray  which  completes  the  reswept  zone  will  probably  be  differ- 
ent than  the  ray  used  and  it  will  have  a  slightly  different 
angle  crossing  the  SLD.   These  angular  errors  will  cause  the 
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greatest  CZ  gain  error  in  the  sin  6~  term  of  the  algorithm. 
Since  sine  is  directly  proportional  to  angle  at  small  angles, 
an  error  of  a  factor  of  two  in  9.  (a  quite  possible  event) 
could  cause  a  gain  error  of  approximately  3  dB . 

Once  CZ  gain  is  computed  and  stored,  the  sound  frequency 
of  interest  is  entered,  and  the  attenuation  coefficient  is 
calculated  using  Thorpe's  equation  (p.  102,  Ref.  1): 


a  =  (0.001094 


2  2 

0.1  f         40  f 


1  +  f2      4100  +  f2 


dB/m)        (6) 


In  Eq.  6,  f  is  in  kHz,  and  the  constant  in  front  of  the  ex- 
pression converts  attenuation  coefficient  from  dB/kyd  to 
dB/m.   The  program  user  enters  frequency  in  Hz,  and  the  pro- 
gram performs  the  conversion  to  kHz. 

Finally,  the  transmission  loss  in  the  n    CZ  annulus 
for  the  frequency  of  interest  (TL  )  is  determined  by: 

TL   =  20  log  (n  RCZi)  +  z  (n  RCZi)  -  G  (7) 

n         3 

In  this  equation,  the  subscript  n  denotes  the  n    CZ  annulus, 

the  range  to  which  is  n  times  RCZi. 

After  a  range  program  is  run,  and  after  the  gain  portion 

of  the  G  and  TL   program  has  been  completed,  TL   values  for 

n    3  n 

a  variety  of  frequencies  and  CZ  annuli  may  be  rapidly  ob- 
tained for  the  SVP ,  source  depth,  and  receiver  depth  condi- 
tions entered.   If,  however,  a  different  set  of  source/receiver 
depth  conditions  are  also  of  interest,  the  entire  procedure 
beginning  with  the  appropriate  range  program  must  be  per- 
formed again. 
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D.   CALCULATOR  PREDICTIONS  COMPARED  TO  ICAPS 

1.   Choosing  the  SVP  Points  for  the  Program 

The  five  point  SVP  limitation  of  the  ray  tracing 
procedure  is  a  rather  serious  handicap  in  many  situations. 
Actual  sound  velocity  profiles  not  only  are  curvilinear  in 
overall  shape  but  also  have  many  small  scale  features  and 
they  are  time  varying  functions.   Approximating  these  curves 
with  only  four  straight  line  segments  presents  a  difficult 
challenge . 

In  general,  matching  the  gradients,  sound  velocities, 
and  associated  depths  are  all  important  in  choosing  SVP 
points.   The  greatest  potential  for  causing  large  range  pre- 
diction errors  occurs  when  the  SVP  contains  an  extensive 
near  surface  layer  with  very  slight  velocity  gradient.   The 
horizontal  distance  traveled  by  a  shallow  depression  angle 
ray  within  such  a  layer  varies  considerably  with  small 
changes  in  the  gradient  or  layer  thickness.   Under  such  con- 
ditions, then,  it  is  extremely  important  to  match  those 
characteristics  as  closely  as  possible. 

Another  important  item  to  carefully  match  is  the 
sound  velocity  at  the  DSC  axis.   This  velocity  determines 
the  maximum  angle  of  depression  for  each  ray  prior  to  com- 
mencing upward  refraction.   The  horizontal  distance  traveled 
by  a  ray  below  the  axis  is  highly  dependent  on  that  angle. 

Table  IV  contains  the  five  point  sound  velocity 
profiles  picked  by  the  author  for  use  in  comparing  the  pro- 
gram performance  to  ICAPS  predictions.   Depths  in  the  table 
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5  Points 

FEB 

MAY 

AUG 

NOV 

PACIFIC 

Dl 

75 

10 

0 

30 

CI 

1489.8 

1500.2 

1508.6 

1497.7 

D2 

340 

80 

100 

100 

C2 

1480.0 

1490.0 

1483.0 

1484.0 

D3 

500 

600 

600 

600 

C3 

1476.2 

1477.0 

1475.5 

1476.0 

D4 

1060 

1600 

1750 

1400 

C4 

1480.0 

1485.0 

1486.5 

1483.0 

D5 

1940 

2610 

3120 

2460 

ATLANTIC 

Dl 

100 

0 

0 

0 

CI 

1524.6 

1530.9 

1541.1 

1535.5 

D2 

550 

100 

125 

125 

C2 

1522.5 

1522.0 

1522.5 

1522.5 

D3 

1080 

575 

550 

625 

C3 

1493.0 

1524.0 

1524.5 

1524.5 

D4 

1750 

1225 

1200 

1175 

C4 

1495.0 

1488.5 

1487.5 

1489.0 

3780 

4175 

4760 

4460 

MEDITERRANEAN 

Dl 

10 

0 

10 

CI 

1526.6 

1537.6 

1528.7 

D2 

30 

50 

50 

C2 

1518.0 

1517.5 

1517.5 

D3 

100 

100 

125 

C3 

1513.0 

1513.0 

1512.0 

D4 

800 

700 

900 

C4 

1521.5 

1519.5 

1522.5 

D5 

1125 

1800 

1290 

Table  IV, 


Five  Point  Sound  Velocity  Profiles. 
(Depths  in  meters,  velocities  in  m/sec) 
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are  in  meters,  and  sound  velocities  are  in  meters  per  second. 
The  reader  may  want  to  plot  these  points  on  the  graphs  of 
Figs.  1(a)  through  1(1)  so  he  may  see  how  the  four  isogradi- 
ent  layers  picked  match  the  ICAPS  profiles.   It  should  be 
pointed  out  that  only  the  initial  selection  of  SVP  points 
was  used  in  the  subsequent  comparisons  of  calculator  model 
results  to  ICAPS  predictions.   Since  an  ASW  aircrewman  using 
the  programs  in  attempting  an  in  situ  prediction  of  acoustic 
conditions  would  not  be  able  to  judge  whether  SVP  point  ad- 
justments would  improve  or  degrade  prediction  accuracy,  it 
was  felt  that  comparing  resuts  of  the  initial  SVP  point  se- 
lection with  ICAPS  would  be  more  meaningful  to  the  objective 
of  developing  the  calculator  programs. 

Comparing  calculator  model  predictions  to  ICAPS  pre- 
dictions in  a  definitive  statistical  manner  was  not  done. 
The  main  reason  for  this  was  alluded  to  in  the  preceding  para- 
graphs.  Since  the  SVP  points  entered  in  the  calculator  pro- 
gram must  be  picked  subjectively  by  the  person  using  the 
program  and  since  it  is  unlikely  different  people  would  pick 
the  exact  same  points  off  any  given  SVP,  it  is  clear  that 
calculator  results  can  be  expected  to  vary  from  operator  to 
operator . 

2 .   CZ  Range  and  Annulus  Width  Comparisons 

The  calculator  range  programs  produce  one  value  each 
for  RCZi  and  RCZo  for  any  given  SVP,  source  depth,  and  re- 
ceiver depth  situation.   Under  the  same  set  of  conditions, 
ICAPS  yields  four  sets  of  RCZi  and  CZW  values,  one  set  for 
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each  of  the  four  frequencies  entered.   In  order  to  compare 
the  calculator  performance  to  ICAPS  it  was  first  necessary 
to  reduce  the  ICAPS  predictions  to  one  value  each  for  RCZi 
and  CZW  for  each  SVP/source/receiver  condition.   This  was 
done  by  simple  averaging  to  eliminate  the  frequency  variable 
from  the  ICAPS  range  and  width  predictions. 

Figures  9(a)  -  9(c)  display  range  and  width  compari- 
sons in  graphical  form  for  the  Pacific,  Atlantic,  and  Medi- 
terranean locations  respectively.  In  each  figure  the  double 
barred  lines  represent  the  ICAPS  first  CZ  annuli  predictions 
(averaged  over  frequency) ,  and  the  single  barred  lines  repre- 
sent the  calculator  predictions.  Numerical  values  for  inner 
and  outer  first  CZ  ranges  may  be  obtained  from  the  scales  at 
the  tops  of  the  figures. 

In  all,  there  were  32  cases  where  these  graphical 
comparisons  could  be  made.   The  following  comments  pertain 
to  those  comparisons: 

a)  In  30  of  the  32  cases  the  calculator  annuli 
overlap  at  least  a  portion  of  the  ICAPS  annuli. 

b)  In  14  of  the  32  cases  the  calculator  annuli 
are  completely  contained  within  the  limits  of  the  ICAPS 
annuli . 

c)  In  all  32  cases  RCZi  ranges  predicted  by 
the  calculator  were  greater  than  those  predicted  by  ICAPS. 
In  the  12  Pacific  cases,  the  calculator  RCZi  values  were  ap- 
proximately 1.7  nm  greater  than  ICAPS  on  the  average.   In 

the  12  Atlantic  cases,  the  average  difference  was  approximately 
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2.8  nm.   And  in  the  eight  Mediterranean  cases,  1.5  nm  was 
the  mean  difference. 

d)   In  27  of  the  32  cases  the  CZ  width  predic- 
tions from  the  calculator  were  more  narrow  than  the  ICAPS 
predicted  widths.   Three  of  the  five  cases  where  calculator 
CZW  exceeded  ICAPS  CZW  were  from  the  February  SVP  in  the  At- 
lantic location.   That  SVP  contained  a  very  deep  (500  meter) , 
nearly  isovelocity  layer  near  the  surface.   In  such  a  pro- 
file, CZ  refraction  produces  ray  paths  that  are  spread  over 
a  very  wide  (in  this  case  16-20  nm)  annulus.   Only  the  rays 
which  return  to  the  SLD  within  the  first  few  nm  at  the  inner 
edge  of  that  annulus  experience  sufficient  convergence  to 
produce  detectable  CZ  gain,  however.   Going  from  inner  to 
outer  edge  of  such  an  annulus  the  CZ  refracted  rays  rapidly 
fan  out  experiencing  progressively  less  convergence  and  pro- 
ducing progressively  less  CZ  gain.   Additionally,  if  the 
bottom  grazing  ray  were  considered,  it  would  be  seen  to  limit 
the  reswept  region  of  this  type  annulus  to  something  far  less 
than  that  indicated.   Since  the  calculator  model  fails  to  ac- 
count for  either  of  these  factors,  it  fails  rather  dramat- 
ically to  produce  a  "practical"  CZ  annular  width  from  this 
SVP  type. 

In  summary,  the  calculator  model  produces  CZ  annuli 
that  roughly  agree  with  those  produced  by  ICAPS  in  all  three 
ocean  basins  considered.   Calculator  RCZi  ranges  are  5-10% 
greater  on  the  average  than  the  ICAPS  values.   Calculator 
CZW  predictions  (excluding  the  Atlantic  February  SVP)  are 
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40-50%  narrower  than  ICAPS  widths  on  the  average.   And,  the 
Atlantic  February  case  indicates  there  is  at  least  one  SVP 
type  in  which  the  calculator  model  fails  to  produce  even 
marginally  acceptable  results  for  CZW. 
3 .   CZ  Gain  Comparisons 

As  with  CZ  range  and  width  comparisons,  it  was  nec- 
essary to  average  the  ICAPS  gain  data  with  respect  to  fre- 
quency before  calculator  gain  predictions  could  be  compared. 
The  estimated  ICAPS  gain  values  in  Tables  II  and  III  were 
thus  reduced  to  one  number  for  each  SVP,  source,  and  receiver 
condition.   Table  V,  CZ  Gain  Prediction  Comparisons,  contains 
numbers  that  represent  the  difference  between  calculator  gain 
predictions  and  the  averaged  ICAPS  values.   Minus  signs  in 
the  table  indicate  those  cases  where  calculator  gain  was  less 
than  the  ICAPS  value. 

As  with  the  range  and  width  comparisons,  the  worst 
agreement  occurred  in  the  Atlantic  winter  SVP  case.   Since 
CZW  is  a  term  in  the  gain  algorithm,  the  extremely  wide  an- 
nuli  predicted  by  the  calculator  caused  gain  values  to  be 
far  too  low  for  the  three  source/receiver  conditions  assoc- 
iated with  that  SVP. 

Excluding  the  Atlantic  winter  SVP  case,  the  follow- 
ing comments  can  be  made  concerning  the  other  29  CZ  gain 
comparisons : 

a)   Calculator  gain  values  ranged  from  7.3  dB 
lower  to  5  dB  higher  than  the  averaged  ICAPS  values. 
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GICAPS    /    ^Gcalc 

GICAPS^ 

Source/ 

SVP 

Receiver 

Profile 

(ft) 

Pacific 

Atlantic 

Mediterranean 

60/60 

13A/-5.4 

15.5/-10.5 

FSB 

60A00 

11  A/-2  A 

13.0/-7.0 

300AOO 

11.2/-2.2 

17.5/-11.5 

60/60 

16. '4/    0.6 

16.5/-0.5 

A. 3/  3.7 

MAY 

60A00 

11 . 6/-2.6 

A.0/-2.0 

10.7/-1 .7 

300AOO 

12.7/-0.7 

15.0/-3.0 

13.3/-6.3 

60/60 

15.8/   2.2 

16.0/  2.0 

17.7/-0.7 

AUG 

60A00 

12.8/-1.8 

11.5/  0.5 

10.0/  0.0 

300/400 

12.4/   0.6 

A.5/-I.5 

60/60 

A.0/-4.0 

12.0/  5.0 

13.0/  4.0 

NOV 

60A00 

11.3/-7.3 

11 .0/   1.0 

11 .0/-3.0 

300A00 

10.7/   1.3 

12.5/   0.5 

13.7/-2.7 

Table  V.   CZ  Gain  Prediction  Comparisons 
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b)  In  22  of  the  29  comparisons,  calculator 
values  were  within  3  dB  of  ICAPS. 

c)  In  nine  of  the  29  comparisons  calculator 
values  were  within  one  dB  of  ICAPS. 

d)  On  the  average,  calculator  gain  values  were 
approximately  one  dB  less  than  ICAPS.   This  result  is  incon- 
sistent with  calculator  CZW  results  in  light  of  the  gain  mod- 
el used.   Since  calculator  CZW  values  averaged  only  slightly 
more  than  half  the  ICAPS  widths,  it  would  have  been  more  con- 
sistent if  calculator  gain  values  turned  out  two  to  three  dB 
higher  than  ICAPS  (acoustic  power  being  spread  over  less  area 
in  the  CZ  annuli,  other  things  being  equal) .   Perhaps  an  ex- 
planation for  this  apparent  discrepancy  is  that  the  calcula- 
tor model  does  not  consider  the  contribution  of  surface 
reflected  energy  adding  to  the  energy  from  upward  traveling 
sound  rays  at  the  receiver  depth.   In  an  actual  CZ  annulus 
the  downward  traveling,  surface  reflected  energy  adds  approx- 
imately three  dB  to  the  CZ  gain  over  much  of  the  annulus 
width.   Apparently,  the  FACT  model  in  the  ICAPS  system  in- 
cludes this  consideration.   It  is  also  apparent  that  neglect- 
ing surface  reflected  energy  in  the  calculator  gain  model  has 
the  effect  of  canceling  errors  that  should  result  from  CZW 
values  being  too  narrow. 

In  summary  of  the  gain  results,  it  can  be  said  that 
the  ray  tracing  technique  used  in  the  calculator  model  worked 
reasonably  well.   Since  three  fourths  of  the  comparison  cases 
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resulted  in  gain  values  within  three  dB  of  the  estimated 
ICAPS  figures,  TL  values  from  the  calculator  displayed  the 
same  close  agreement. 
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IV.   CONCLUSIONS 

A.   LIMITATIONS  OF  THE  MODEL  DEVELOPED 

The  HP-67/9  7  calculators  used  in  programming  the  CZ  pre- 
diction model  were  stretched  to  their  limits  in  both  data 
storage  and  program  step  capacity.   Although  not  known  for 
certain,  the  author  feels  significantly  more  accurate  results 
would  be  possible  from  a  calculator  with  only  moderately 
larger  storage  capacity. 

The  data  storage  limitation  which  allowed  only  five  SVP 
points  to  be  entered  is  quite  restrictive  and  no  doubt  plays 
a  large  role  in  the  CZ  range  and  width  inaccuracies  obtained. 

Program  step  capacity  forced  several  short  cuts  to  be 
taken  which  again  would  not  have  been  necessary  with  a  mod- 
erately larger  program  memory.   Two  separate  range  and  width 
programs  were  required  due  to  insufficient  program  space  to 
incorporate  tests  for  different  source  and  receiver  depth 
cases.   Also,  source  and  receiver  depths  are  strictly  allowed 
only  within  the  upper  two  SVP  isogradient  layers  because  pro- 
gram space  was  not  available  to  check  for  the  correct  layer 
if  all  depths  were  allowed.   Additionally,  and  perhaps  the 
greatest  source  of  CZW  errors  observed,  program  step  limita- 
tion prevented  incorporating  a  method  of  considering  the 
bottom  limited  CZ  sound  ray  in  determining  the  range  to  the 
outer  edge  of  a  CZ  annulus .   The  program  developed  ignores 
the  bottom  entirely  and  considers  only  the  reswept  zone  in 


predicting  annular  width.   Since  calculator  CZW  results  were 
considerably  shorter  than  those  indicated  by  ICAPS ,  it  is 
assumed  the  discrepancy  is  due  to  not  considering  CZ  rays 
beyond  the  reswept  zone.   The  first  priority  in  making  im- 
provements to  the  calculator  model,  should  a  larger  capacity 
machine  be  implemented,  would  be  incorporating  a  better  meth- 
od for  selecting  the  ray  which  defines  the  outer  limit  of 
the  CZ  annulus . 

B.   USEFULNESS  OF  THE  MODEL  DEVELOPED 

The  degree  of  success  in  producing  a  useable  CZ  prediction 
model  for  handheld  calculators  must  be  determined  by  consid- 
ering the  objectives  set  forth  in  the  first  section  of  this 
study.   The  central  idea  was  to  ascertain  if  a  calculator 
model  would  improve  on  ASRAPS  CZ  prediction  accuracy  in  the 
case  where  3T  conditions  determined  in  situ  differed  from 
those  used  to  generate  the  ASRAPS  TL  profiles.   Inherent  in 
this  objective  is  the  assumption  that  ASRAPS  TL  profiles 
generated  primarily  from  climatological  data  would  be  in 
error  due  to  lack  of  input  data  accuracy.   Also  inherently 
assumed  is  that  given  identical  input  data  the  calculator 
model  would  produce  less  accurate  results  than  the  digital 
computer  model(FACT)  used  in  ASRAPS(and  ICAPS)  due  to  obvious 
differences  in  data  and  program  capacities.   The  real  question 
then  is  a  trade-off  comparison:   Will  the  basically  less  accu- 
rate calculator  model  produce  better  CZ  predictions  with  actual 
environmental  data  than  the  more  sophisticated  digital  com- 
puter model  which  had  only  climatological  input  data? 
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Before  addressing  the  answer  to  this  question,  character- 
istics of  the  calculator  model  developed  will  be  compared  to 
the  list  of  six  desirable  characteristics  described  in  sec- 
tion I . 

1.  Easily  available  input  data. 

The  data  required  are  an  SVP,  assumed  source 
depth,  hydrophone  depth,  and  frequency  of  interest.   The 
only  portion  of  this  information  not  presently  available  to 
ASW  aircrews  is  that  part  of  the  SVP  below  the  1,000  ft 
depth  limit  of  the  AN/SSQ-36  bathythermograph  buoy.   SVP  data 
from  the  surface  to  1,000  ft  (the  area  where  seasonal  and 
diurnal  variations  predominantly  occur)  is  easily  obtained 
from  the  BT  buoy  information. 

2.  Ease  of  program  operation. 

Anyone  familiar  with  HP-67/9  7  calculator  use 
could  operate  this  program  without  additional  training. 

3.  Output  data. 

The  program  provides  CZ  annulus  range  and  width 
as  well  as  TL  values  for  all  frequencies  of  interest  in  all 
annuli  of  interest. 

4.  Short  run  time. 

To  run  a  complete  program  requires  approximately 
10  minutes  once  SVP  data  is  obtained.   Deploying  a  BT  buoy 
and  converting  the  temperature  trace  to  an  SVP  would  take  an 
additional  10-15  minutes. 
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5.  Based  entirely  on  acoustic  theory. 

The  program  uses  only  ray  tracing  techniques  in 
producing  its  output  terms. 

6 .  Agreement  with  large  computer  models . 
Calculator  CZ  ranges  obtained  averaged  5-10% 

greater  than  ranges  obtained  from  ICAPS.   CZ  width  results 
averaged  only  40-50%  of  those  obtained  from  ICAPS.   And, 
there  was  one  SVP  case  studied  (winter,  Atlantic)  in  which 
the  calculator  CZW  results  were  very  different  from  ICAPS. 
That  SVP  case  was  considered  a  failure  of  the  calculator 
model,  and  it  must  be  conceded  the  model  does  not  work  for 
all  CZ  situations.   Excluding  the  obvious  CZW  failure  SVP 
case,  TL  values  from  the  calculator  averaged  about  one  dB 
lower  than  ICAPS  with  extreme  deviations  observed  ranging 
from  -7.3dB  to  +5dB  around  the  ICAPS  values.   Additionally, 
calculator  results  can  be  expected  to  vary  from  operator  to 
operator  since  SVP  points  must  be  picked  subjectively  from 
an  SVP  graph. 

Returning  to  the  main  objective  of  the  study,  the  author 
feels  that  only  half  of  the  trade-off  question  has  been 
answered.   An  easily  operated,  purely  theoretical  model  was 
developed  which  works  for  most  (but  not  all)  CZ  producing 
SVP  conditions.   And  a  measure  of  its  accuracy  compared  to 
the  sophisticated  FACT  computer  model  was  obtained.   Yet  to 
be  answered  is  how  inaccurate  ASRAPS  CZ  predictions  are  when 
observed  BT  conditions  differ  from  climatological  conditions 
This  portion  of  the  question  is  very  difficult  to  answer  and 
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in  fact  would  be  a  very  large  study  in  itself.   Generally, 
inaccuracies  must  range  over  a  scale  from  insignificant  to 
considerable  as  environmental  deviations  range  from  slight 
to  great.   The  most  likely  variables  affecting  degree  of  in- 
accuracy are  surface  water  temperature,  mixed  layer  depth, 
and  in  layer  and  below  layer  gradients.   The  effects  of  vary- 
ing these  or  other  possible  factors  one  at  a  time  or  in 
various  combinations  on  CZ  range,  width,  and  gain  must  be 
known  before  the  entire  question  can  be  answered.   Further, 
definite  magnitudes  of  environmental  factor  deviation  must 
be  determined  so  that  a  person  can  judge  when  ASRAPS  inaccu- 
racies are  likely  to  be  greater  than  the  calculator  model 
inaccuracies.   Until  these  points  are  answered  it  would  be 
inappropriate  to  recommend  use  of  the  calculator  model  as  a 
routine  method  of  updating  ASRAPS  CZ  predictions  in  situ. 
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APPENDIX 

HP-67/9  7  Calculator  Programs  for  Convergence 
Zone  Range ,  Width  and  Transmission  Loss  Predictions 


Steps  required  to  use  the  programs: 

1.  Deploy  a  bathythermography  buoy  in  the  operating  area 
of  interest. 

2.  Convert  the  BT  buoy  information  to  a  sound  velocity 
profile  of  the  upper  1,000  ft  of  the  ocean  area. 

3.  Combine  the  upper  SVP  data  with  a  graph  of  climato- 
logical  SVP  data  which  depicts  sound  velocity  conditions 
below  the  1,000  ft  level. 

4.  Pick  five  points  from  the  combined  SVP  graph  which 
best  represent  the  deep  sound  channel  portion  of  the  SVP. 

The  first  point  should  be  at  the  sonic  layer  depth,  the  fifth 
point  at  the  bottom  of  the  DSC  where  sound  velocity  equals 
that  at  the  SLD,  and  the  other  three  points  at  points  on  the^ 
graph  such  that  when  straight  lines  are  drawn  to  connect  the 
five  points  they  create  a  linearly  segmented  SVP  which  matches 
the  actual  SVP  as  closely  as  possible. 

5.  Pick  the  appropriate  CZ  Range  and  Width  program  to  be 
used  as  follows:   If  both  source  and  receiver  are  below  the 
SLD  use  the  Deep/Deep  program.   If  source  or  receiver  or  both 
are  above  the  SLD  use  the  Shal/Shal  or  Crosslayer  program. 

6.  Load  and  run  the  appropriate  Range  and  Width  program 
according  to  the  accompanying  instructions. 

7.  Leaving  the  calculator  power  on  and  data  storage  reg- 
isters unchanged,  load  and  run  the  Gain  and  Transmission  Loss 
program  according  to  its  instructions. 


A  word  about  units : 

As  currently  written,  the  programs  use  metric  units; 
meters  for  depths,  and  m/sec  for  sound  velocities.   To  con- 
vert the  programs  for  english  unit  input  data,  feet  for 
depths  and  ft/sec  for  velocities,  the  conversion  factor 
1,852  m/nm  should  be  changed  to  6,075  ft/nm  where  occurring. 
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User  Instructions  for  CZ  Range  Programs 
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Input 
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Output 
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liminary calculations 

Enter  Receiver  Depth  and 
Source  Depth. 

Range  to  inner  edge  of 
CZ  (RCZi)  and  range  to 
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of  User  Control  keys  B 
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Storage  Allocation  for  CZ  Range  Programs 
Registers : 


RO 
Rl 
R2 
R3 
R4 
R5 
R6 
R7 
R8 
R9 


D   /  ZAr 


D2  /  CS 


D3  /  CR 


D.  /  r  .   /  RCZi 

4  '       mm  ' 


D5  /  rQ  /  RCZo 


'0 


SO 
SI 
S2 
S3 
S4 
S5 
S6 
S7 
S8 
S9 


Initial  Flag  Status  and  Use 


^2 
^3 
^4 


0:   Off,  Unused 


Off,  On  if  RCVR  is  shallow 
in  Shal/Shal  or  Crosslayer 
program  only. 


S: 

C^  /  cos  e 

0        rmin 

B: 

A      /  A 

SLD  '       rmin 

C: 

cos  9 
w 

D: 

cos    .  , 

0W+1 

E: 

r 

w 

I: 

Control 

Off,  set  prior  to  r 


0 


calculation  in  Deep/Deep 
program  only. 

Off,  set  by  data  entry 
until  r_  is  found  in 

S/S  or  Crosslayer  program 
only. 


Display  Status:   DSP  1 


User  Control  Keys 


Data  entry 
Display  RCZi 
Display  RCZo 


a 
b 
c 
d 
e 
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Step  Key 

s   Code 

Explanation 

Ste 

P   Key 

s   Code 

Explanation 

003 

♦LdLh 

570i' 
ISZI 

21  11 

35  45 

It)  2b   4b* 

Enter  data 

05J 

f?/S 

51 

Enter  DR f  DS 

052 
053 

RCLB 

36  68 
-45 

m 

R/S 

51 

854 
055 
056 

FCL  ', 

x 

RCLl 

36  45 

-35 

36  01 

Compute  Co 

O 

885 

000 

0 

-62 

00 

Compute 
Gradients 

00  7 

£ 

02 

s0 

057 

+ 

-55 

00S 
00? 

rszi 

35  45 
16  26  46 

058 

ST  02 

35  8c 

055 

060 

Ri 

RCL8 

-31 
36  88 

010 

RCL  3 

36  03 

tfli 

012 

RCLl 

RCL2 

36  01 
3b  02 

gl 

861 
862 

RCL  i 

-45 
36  45 

Compute  C0 

813 

RCL0 

36  00 

063 

X 

-35 

014 

QSBa 

23  16  Jl 

064 
065 

RCLl 

+ 

36  01 
-55 

"015 

RCL5 

36  05 

016 
017 

010 
019 

RCL3 
RCL4 
RCL2 
QSBa 

36  03 
36  04 
36  02 

23  16  11 

g2 

066 

ST04 

35  84 

067 

866 
069 

SF2 

1 

S70C 

16  21   62 
61 

35  13 

Initialize 
Ar  routine 

020 

RCL7 

36  07 

878 

*LBL2 

21  82 

021 

RCL5 

36  05 

071 

asBc 

23  16  13 

022 

RCL6 

36  06 

s3 

072 

RCL3 

36  83 

023 

RCL4 

36  04 

073 

PCLC 

36  13 

2  A*! 

024 

QSBa 

23  16  11 

074 

075 

QSBai 
ST  08 

23  16  14 
35  00 

025 

RCL1 

36  01 

02b 

027 

RCL7 
RCL8 

36  07 
36  08 

% 

076 

ST+8 

35-55  00 

077 

RCL3 

36  83 

028 

PCL6 

36  06 

078 

EHTt 

-21 

023 

QSBa 

23  16  11 

879 
080 
881 

EHTt 
RCL5 
RCLD 

-21 
36  85 
36  14 

2  Ar2 

039 

eset 

23  16  12 

Set   1=10 

831 

RCL0 

36  00 

032 

X*0? 

16-42 

Mixed  layer  ? 

882 

QSBal 

23  16  14 

033 

ewe 

22  08 

083 
084 

ST+8 

$7+8 

35-55  00 
35-55  00 

034 
035 

ST  08 
QT01 

35  08 
22  01 

If  no  ML, 
0   to   R8 

085 
086 

RCL  5 

Em-r 

36  05 

-21 

836 

*LBL8 

21  80 

037 

RCL9 

36  89 

887 

EHTt 

-21 

030 

X 

-35 

888 

RCL7 

36  87 

2  A^? 

039 

CHS 

.  -22 

co 

089 

RCLD 

36  14 

040 

RCL1 

36  81 

898 

QSBd 

23  16  14 

041 

+ 

-55 

891 

ST+8 

35-55  80 

042 

STOrt 

35  11 

892 

ST+8 

35-55  00 

043 

esse 

23  16  13 

893 

RCLl 

36  81 

044 

PCLh 

36  11 

894 

RCLD 

36  14 

045 

1 

01 

ML  present, 

895 

T 

-24 

04c 

es&orf 

23  16  14 

2  Ar0    to   R8 

896 

RCL ; 

36  45 

0       a  r> 

047 

ST08 

35  08 

897 

T 

-24 

2  Ar>4 

84S 

ST+8 

35-55  08 

898 
899 

100 

RCLD 
GSffe 

36  14 

23  16  15 

-35 

049 

050 

tLBLl 
GSBk 

21  81 
23  16  12 

Set   1=10 

Z   Range  Program    (Deep/peep   Case) 
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Step  Key 

3   Code 

Explanation 

Ste 

P    Key 

s   Code 

Explanation 

J0J 

182 

ST+8 

sr+e 

35-55  88 
35-55  88 

151 
152 

153 

9SBd 
ST+8 

DSZI 

23 
16 

16    14 

-55  88 

25  46 

to   rQ 

163 

GSBit 

23  16  12 

Set   1=10 

184 
185 
IBS 

F2? 
GTQ3 

QT04 

16  23  82 
22  83 
22  84 

rn   found   ? 

154 
155 
156 
157 
15S 

QSBc 
RCL4 
RCL* 

GSBd 
ST-6 

23 
35- 

16  13 
36  84 
36  11 
16  14 
-45  86 

to   r    . 
mm 

197 

iee 

109 

lie 
in 

112 

*LBL3 
RCL8 
ST+8 

1 

+ 

STOE 

21  63 

36  88 

35-55  88 

81 

-55 

35  15 

Store   rQ 

&   first   r 
w 

159 
168 
161 
162 
163 
164 

1 
8 
5 

STrQ 

35- 
35- 

61 
68 
05 
82 
-24  86 
-24  88 

Convert 

RCZi   &  RCZo 
to  nm 

113 

114 

115 
116 
117 

*LBL4 
RCL8 
RCLE 

QT05 

21  84 
36  88 
36  15 
1 6-35 

22  85 

r    .      found   ? 
mm 

165 

essfc 

23 

16  12 

Set   1=10 

166 

167 
16B 

*LBLB 

RCL6 

R/S 

21   12 
36  86 

51 

Display  RCZi 

118 

119 

Ri 
STOE 

-31 
35  15 

Store 
next   rw 

169 
176 

171 

tLBLC 

RCLB 

R/'S 

21   13 

36  88 

51 

Display  RCZo 

128 
121 
122 
123 
124 
125 
126 

Kft.fi 

5 

+ 

ST  OB 

COS 

STOC 

36  12 

-62 

85 

-55 

35  12 

42 

35  13 

Increment 

9SLD  & 
cos   9SLD 

172 
173 
174 
175 
176 
$177 
178 
179 
136 

*LFLa 

Ri 

Rt 

STQi 
1SZI 

P7,S 

21 
16 

16  11 
-45 
-31 
-45 

16-31 
-24 

35  45 

26  46 
24 

Gradients 
Subroutine 

121 

QT02 

22  82 

Next   r  rout in 

J26 
129 

*LBL5 
ST06 

21  85 
35  86 

Store   r    . 

mm 

130 
131 
132 
133 
134 
135 

RCLB 

• 

5 

COS 
STOfi 

36  12 
-62 

85 

-45 

42 

35  11 

Initiate 
corrections 

181 
182 
183 
184 
185 

tLBLh 

1 

8 

ST  01 

RTN 

21 

16  12 
61 
86 

35  46 

24 

Set   1=10 
Subroutine 

136 
137 
13$ 
139 
148 
141 

QSBc 
RCL2 

1 
GS8J 
ST+8 
OSZI 

23  16  13 
36  82 
81 
23  16  14 
35-55  88 
16  25  46 

corrections 

to   rQ 

186 
187 
188 
189 
198 

*LBLc 

RCL1 

EtiTt 

EHTt 

RTH 

21 

16  13 
36  81 

-21 

-21 
24 

Initiation 
Subroutine 

191 
192 
193 
194 
195 
196 
197 
198 
199 

26e 

*LEU 
STOC 

X 

ST  00 
GSBe 
RCLC 

C-SBe 

21 

16  14 

35  13 

-35 

-41 

-24 

35  14 
16  15 

36  13 
16  15 

-45 

A^ 

Subroutine 

142 
143 
144 
145 
146 
147 

bSBc 
RCL2 
RCLh 
*$BJ 
ST -6 
DSZI 

23  16  13 
36  82 
36  11 
23  16  14 
35-45  86 
16  25  46 

to   r 

mm 

148 
149 
158 

QSBc 

RCL4 

i 

23  16  13 

36  84 

61 

corrections 

CZ  Range  Program  (Deep/Deep  Case) 
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Step  Keys   Code 


281 
282 
222 
264 
225 
226 
227 
288 


289 
218 
211 
212 

213 
214 
215 
216 


x 
RCLC 

PCLi 
ISZI 

RTH 


-35 

36  13 

-24 

36  45 

16  26  46 

-24 

16  31 

24 


tLSLe 

CHS 
1 

+ 

RTH 
R/.S 


21  16  15 
53 

-22 
81 

-55 
54 
24 
51 


Explanation 


A** 

Subroutine 

(continued) 


cos  to  sin 
Subroutine 


J tep  K  ey s  J  ode 


Explanation 


CZ  Range  Program  (Deep/Deep  Case) 
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Step  Key 

s   Code 

Explanation 

St€ 

ip    Ke.Y 

s   Code 

Explana tion 

B81 

m 

883 

884 

*l£Lh 
iTQi 
JSZ1 

16 

21   11 

35  45 

26  46 

51 

inter  data 

05  i 
852 
857 
854 
855 

RCL9 
RCLl 
8702 

Jo  89 

-35 

36  81 

-  JlJ 

35  82 



885 

886 

887 
m 

88S 

8 
c 

STOi 

JSZI 

16 

-62 
88 

82 
35  45 
26  46 

Compute 
Gradients 

§0 

856 
857 
858 
859 
868 
861 
862 
863 
864 
865 
866 
867 
866 
869 
878 

Ri 

RCl8 

X>  V""' 

QT08 

GT01 

*LBL8 

SF1 

DSZI 

*LBL1 

RCL  i 

x 

RCL  J 
+ 
ST04 

36  88 
16-34 
22  88 
22  81 
21  88 
16  21  81 
16  25  46 
21   81 

-45 
36  45 

-35 
36  81 

-55 
35  84 

Set  Flag  0 
for  shallow 
Receiver 

and 

Compute  C^ 

010 

en 

812 
813 

814 

RCL3 
RCL1 
RCL2 
RCL8 

C-SBa 

23 

36  83 
36  81 
36  82 
36  88 

16  11 

% 

815 
81€ 
817 
81S 

813 

RH5 
RCL3 
RCL4 
RCL2 
GSB<x 

23 

36  85 
36  83 

36  84 
36  82 

16  11 

°2 

82* 
821 
822 
823 
824 

RCL7 
PCL5 
RCL6 
RCL4 

bSBa 

23 

36  87 
36  85 
36  86 
36  84 

16  11 

g<3 

871 

GS66 

23  16  12 

Set   1=10 

872 
873 
874 
875 
876 
877 
879 

*LBL2 

QSBc 
RCL3 
FXLC 
QSBoi 
£708 
ST+8 

21   82 
23  16  13 

36  63 

36  13 
23  16  14 

35  88 
35-55  88 

2  A^1 

825 
826 

827 
828 
829 

RCL1 
RCL7 
RCL8 
RCL6 
GSBa 

23 

36  81 
36  87 
36  86 
36  86 

16  11 

^ 

879 
888 
8S1 
882 
883 
884 
885 
886 

RCL3 

EtiTi 

enti 

RCL5 
RCLD 
GSBd 

ST+8 

ST+8 

36  83 

-21 

-21 

36  05 

36  14 

23  16  14 

35-55  88 

35-55  88 

2  Ar2 

838 
831 

STOI 

89 
35  46 

Set   1=9 

832 
833 
834 
835 
836 
837 
838 

RCL8 
RCL9 

X 

CHS 
RCLl 

+ 
ST08 

36  88 

36  89 

-35 

-22 
36  81 

-55 
35  11 

Compute  C 

887 
888 
889 
898 
891 
892 
893 
894 

RCL5 
EHTt 
EHTt 
RCL7 
RCLD 

esBd 

ST+8 
ST+8 

36  85 

-21 

-21 

36  87 

36  14 

23  16  14 

35-55  00 

55-55  00 

2  A^3 

839 
848 
841 
842 
843 
844 

GSEc 
RCLh 
1 
QSBd 
ST08 
ST+8 

23 

23 
35 

16  13 
36  11 
81 
16  14 
55  88 
-55  96 

2  ^r      to   R8 

845 

esBb 

23 

16  12 

Set   1=10 

895 
896 
897 
898 
899 
188 

RCLl 
RCLD 

RCL  i 

RCLD 

56  81 

36  14 
-24 

36  45 
-24 

36  14 

2   A^ 

"  846 

R/S 

51 

inter  DR      DS 

847 
848 

X>  Y? 

k:y 

16-34 

-41 

Reciprocity 
test 

849 
858 

RCL8 

36  88 
-45 

Compute   C„, 

Z  Range  Program  (Shal/Shal  &  Crosslayer  Cases) 
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Step   Key 

s   Code 

Explanation 

Step    Ke.v 

S     1 

:ode 

Explanation 

101 

182 

193 
194 

QSBe 

X 

ST+8 

sr+a 

23  16  15 

-35 

35-55  88 

35-55  88 

151 
152 

154 
155 
1 56 
157 
158 
159 

DSZI 
QSBc 
RCL4 
I 
QSBd 
DSZI 

Fl  ? 

CHS 
ST+8 

16 
23 

23 
16 
16 

35- 

25  46 
16  13 
36  84 
81 
16  14 
25  46 
23  81 

-55  88 

corrections 

to   rQ 

195 

GSSk 

23  16  12 

Set   1=10 

196 
197 
188 

F3? 
ST03 
GT04 

16  23  83 
22  83 

22  84 

rQ   found   ? 

189 

119 
111 
112 

113 
114 

fLSL3 
RCL8 
ST+8 

1 

+ 

STOE 

21  83 

36  88 

35-55  88 

81 

-55 
35  J  5 

Store   rQ 

&  first   r 
w 

169 
161 
162 
163 
164 
165 
166 

QSBc 
RCL4 
RCL8 
QSBd 

Fl? 

CHS 
ST-6 

23 
16 

35- 

16  13 
36  84 
36  11 
16  14 
23  81 
-22 
-45  86 

to   r   . 

mm 

115 
US 
117 
118 
119 

*LBL4 
RCL8 
RCLE 
X&Y? 
QT05 

21  84 
36  88 
36  15 
16-35 

22  85 

r    .      found   ? 
mm 

167 
168 
169 
178 
171 
1  72 

1 
8 

5 

£ 

ST  to 

STr8 

35 
35- 

81 
88 
85 
82 
-24  86 
-24  88 

Convert 

RCZi   &  RCZo 

to  nm 

129 
121 

Ri 
STOE 

-31 
35  15 

Store 
next   rw 

122 
123 
124 
125 
126 
127 
128 

RCLE 

5 

+ 

ST  OB 

COS 

STOC 

36  12 

-62 

85 

-55 

35  12 

42 

35  13 

Increment 

9SLD  & 
cos   G3LD 

1 73 

GSBb 

16  12 

Set   1=10 

174 

175 
176 

*LEL6 

RCL6 

R'S 

21   12 

36  86 

51 

Display  RCZi 

1  »   r 

178 

ei79 

*LBLC 

RCLS 

R  S 

21   13 

36  88 

51 

Display  RCZo 

129 

5702 

22  82 

Next   r   rout in 

138 
131 

*LBL5 

ST  06 

21  85 
35  86 

Store   r    . 
mm 

188 
181 

1S2 
183 
184 
185 
186 
187 
IBS 

tLBLa 
Ri 

Rt 

sto: 

ISZl 
RTH 

21 
16 

i6  11 
-45 
-31 
-45 

16-31 
-24 

35  45 

26  46" 
24 

Gradients 
Subroutine 

132 
133 
134 
135 
136 
137 

RCLB 

5 

COS 
STOk 

36  12 

-62 

85 

-45 

42 

35  11 

Initiate 
corrections 

'  138 
139 
148 
141 
142 
143 
144 

QSBc 

RCL2 

I 

DSZI 
QSBd 
DSZI 

ST+8 

23  16  13 
36  82 
81 
16  25  46 
23  16  14 
16  25  46 
35-55  88 

corrections 
to   rQ 

189 
198 
191 
192 
193 

*LBLh 

1 

8 

STO  I 

RTH 

21 

Jt>  12 
81 
89 

35  46 

24 

Set    1=10 
Subroutine 

194 
195 
196 

197 
198 

tLBLc 

RCL1 
EtiTt 

EHTt 
RTH 

21 

16  13 
36  81 

-21 

-21 

24 

Initiation 
Subroutine 

145 
146 
147 
148 
149 

b SBc 
RCL2 
RCLh 
QSBd 

ST +6 

23  16  13 

36  82 

36  11 

23  16  14 

35-55  86 

to   r    . 
mm 

199 
288 

xLBLj 

STOC 

21 

16  14 
35  13 

158 

Fl? 

16  23  81 

CZ  Range  Program  (Shal/Shal  8c   Crosslayer  Cases) 
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Step 


Code 


29 1 
292 
263 
294 
295 
296 
29? 
296 
299 
219 
211 
212 
213 
214 
215 
216 


217 
218 
219 
229 
221 
222 
223 
224 


x 
X2Y 

ST  00 
QSBe 
RCLC 

QSBe 

x 
RCLC 

RCLi 
I  SI  I 

RBS 

RTH 


23 


-35 

-41 
-24 

35  14 
16  15 

36  13 
23  16  15 

-45 

-35 

36  13 

-24 

36  45 

16  26  46 

-24 

16  31 

24 


tLBLe 

X* 

CHS 

1 

+ 

■fX 
RTH 

g/S 


21 


16  15 

53 

-22 

91 

-55 

54 

24 

51 


Kxplana  tion 


Subroutine 


cos  to  sin 
Subroutine 


tep    Keys    J  ode 


;Xpla  r  i  a  t,  1  ori 


CZ  Range  Program  (Shal/ohal  &  Grosslayer  Cases) 
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User  Instructions  for  CZ  Gain  and  Transmission  Loss  Program 


Step 


Instructions 


Input    Keys 


Output 


Run  one  of  the  CZ  Range 
Programs .   Leave  the 
calculator  on  and  all  of 
the  data  storage  regis- 
ters unchaged. 

Load  the  CZ  Gain  and 
Transmission  Loss  Program. 


Enter  a  zero  if  the  Shal/ 
Shal  or  Crosslayer  program 
was  used,  or  enter  a  one 
if  the  Deep/Deep  program 
was  used. 


0 

or 

1 


G  (dB) 


Enter  a  frequency  of  in- 
terest in  Hz. 


f  (Hz)  B 


a  (dB/m) 


Enter  the  number  of  a  CZ 
annulus  of  interest .  (Step 
5  may  be  repeated  for  as 
many  CZ  annuli  as  desired.) 


n 


TL   (dB 
n 


To  determine  the  TL 


n 


values  for  different 
frequencies ,  return  to 
step  4 . 
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Storage  Allocation  for  CZ  Gain  and  Transmission  Loss  Program 


Reg 

isters : 

RO: 

Ar 

SO: 

gl 

A: 

9S0  /  9R0  /  f2 

Rl: 

ci 

SI: 

g2 

B: 

9SLD  /  9rswp 

R2: 

cs 

S2: 

g3 

C: 

cos  8   / 
w 

R3: 

C2 

S3: 

g4 

cos  9 

rswp 

R4: 

CR 

S4: 

a  terms 

D: 

cos  9  , , 
w+1 

R5: 

C3 

S5: 

E: 

r0  /  9SR  /  9RR 

R6: 

RCZi 

S6: 

I: 

Control 

R7: 

C4  / 

a 

S7: 

R8: 

RCZO 

S8: 

R9: 

G  terms 

/  G 

S9: 

Initial  Flag  Status  and  Use 

0:   Off,  On  for  shallow 
sound  source. 


2:   Off,  unused 


Off,  On  to  decrease  the 

increments  used  in 

finding  9^T  ^  . 
3   SLD 


3:   Off,  set  by  data  entry 
until  rn  found. 


Display  Status:   DSP  0 


User  Control  Keys: 


A: 

Compute  G 

B: 

Compute  a 

C: 

Compute  TL 

D: 

E: 

n 


a: 

b: 
c : 
d: 
e : 
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Step   Key 

s   Code 

Explanation 

Step    Iveyn   Code 

Sx plana  I  Lot; 

00  J 

862 

003 

SF0 

21  ii 

16-43 

16  21   00 

Set   FO    for 
Shal   Source 

051 
052 

053 

1 

+ 

QT03 

61 

-55 

22  63 

Lar§e    9SLD 
increment 

00.* 
005 

1 
STOC 

01 

35  13 

Initialize 
r  routine 

054 
055 
056 
057 
056 
055 
060 

tLSLl 

SFi 
RCLB 

Q 
QTQ3 

21  61 
16  21   01 

36  12 

-62 

05 

-45 

22  83 

First    small 

increment 

of   0„Tn 

906 

00? 
00? 
010 

en 

612 

8l3 

W.SL0 

FCL1 
ENTf 
£N7f 

RCL3 
RCLC 
GSBi 
ST  06 

21      00 

36  01 
-21 
-21 

36  03 

36  13 
23  16  11 

35  00 

2  Arl 

661 
662 

663 

*LBL2 

x>y? 

&TC4 

21  62 
16-34 

22  04 

©             found   ? 
rswp 

014 
815 
816 
617 

010 

015 
020 

RCL3 
EHTt 
EHTt 
RCL5 
RCLD 
QSBa 
ST+8 

36  03 
-2i 

-21 

36  05 

36  14 

23  16  11 

35-55  00 

2  A?o 

664 
065 
666 
067 

RCLB 

1 

36  12 

-62 

91 

-55 

Small   &-,T^ 
increment 

068 
065 
070 

971 
972 

$LBL3 

ST  OB 

COS 

STOC 

GTO0 

21  03 
35  12 

42 
35  13 

22  66 

Store   new 

9SLD   & 

COS       ©oTT, 

021 
022 
023 
024 
025 
026 
027 

RCL5 
EHTt 
EHTt 
RCL7 
RCLD 

ese<x 

ST+8 

36  05 

-21 

-21 

36  07 

3b   14 

23  16  11 

35-55  66 

2  A^3 

973 

674 
075 
676 

677 
67 8 
679 

$L3L4 
1 
S 
5 

ST:<6 

ST*8 

21   64 
61 
06 
05 
02 
35-35  06 
35-35  08 

Convert 

RCZi   3c   RCZo 
to  meters 

028 
029 
030 
031 
032 
033 
634 
035 
036 
037 

RCL1 
RCLD 

RCLS 

RCLD 
QSBb 

X 

ST+6 
ST+6 

36  01 
36  14 

-24 
36  45 

-24 
36  14 

23  16  12 

-35 

35-55  00 

35-55  00 

2   4^ 

686 
631 
682 

RCLB 

COS 

STOC 

36  12 

42 

35  13 

Store 

cos   G__ 

rswp 

683 
684 
885 
986 
687 
883 

RCL2 

X 

RCL1 

COS-' 
STOE 

36  02 
-35 

36  01 
-24 

16  42 

35  15 

&SR 

033 
639 
040 

i 

6 

ST  01 

01 
00 

35  46 

Set    1=10 

889 
090 
891 
892 
893 
894 
895 

RCL2 
RCL1 

cos- 
Fd* 
CHS 

ST  Oh 

36  62 

36  01 

-24 

16  42 

16  23  00 

35  11 

6S0 

(negative    if 
source    shal) 

041 
042 
043 

RCL6 
F3? 

STOE 

36  06 

16  23  03 
35  15 

Store   rQ 

644 

645 
646 

647 

RCLE 

Fl? 

£702 

36  15 

-41 

16  23  61 

22  62 

Check   for 
Small    inc . 

of   93LD 

896 
897 
898 

0-R 

ST  09 

-55 
16  45 
35  89 

A© 

646 
649 

X)  V'V 
GT01 

1 6-34 

22  61 

Decrease 
©3LD   inc.    ? 

899 
188 

RCLE 

RlLh 

36  15 
36  11 

050 

PCLB 

36  12 

)Z  Cain  and  Transmission  Loss  Program 
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Step   Keys   Code 

Explanation 

Step   Keys  Code 

Sxplana  t, ion 

J0j          -               -45 

182  2               82 

183  t                -24 

184  COS                42 

185  ST*9     35-35  89 

cos   ©1 

151  4                 64 

152  1                61 

1 53  0                00 

1 54  0                   00 

155  RCLR          36  11 

156  +                -55 

157  r                   -24 
J  58          1+                56 

Second   term 

186  RCLC          36  13 

187  RCL4           36  84 
1 08          x                -35 
189      PCLl           36  81 
118            r                   -24 
111       COS"          16  42 

9RR 

159      RCL2          36  56 

166       ST 07          35  67 

Store  a 

161           2-           16  56 

(Clear  S4) 

112  RCL4          36  84 

113  RCLl           36  61 

114  t                -24 

115  COS-'          16  42 

9R0 

162  1                 8l 

1 63  .               -62 

1 64  0                    00 

1 65  9                69 

1 66  4                04 

167  EEX              -23 

1 68  3                63 

1 69  CHS              -22 
176       STx7     35-35  87 

Convert  a 

from   dB/kyd 

to   dB/m 

116  +                -55 

117  2                82 

US           r                 -24 

G2 

119         SIN                41 
128      STr9     35-24  69 

sin  9 2 

121  RCLS          36  83 

122  RCL6          36  86 

123  STx?     35-35  89 

RCZi 

17 1  RCL7           36  67 

172  R  S                51 

Display   a 

173  *LBLC          21   13 

174  RCLS          36  66 

1 75  x                -35 

176  EH1T               -21 

nRCZi 

124  -45 

125  ST±9     35-24  89 

CZW 

126  RCL9          36  83 

127  LOG           16  32 

128  1                81 

129  8                86 
136          x                -35 

131  ST 09          35  69 

132  R/S                51 

Compute  and 
Display  G 

177  lOQ           16  32 

1 78  2                82 

1 79  8                86 
1 36          *                 -35 

201og(nRCZi) 

lbl         ;<:'f               -41 
1S2      RCL7          36  87 
183                           -35 

_a  (nRC  Z_i )  _ 

133  *LBLB          21   12 

134  EEX              -23 

135  3                83 

136  t                -24 

Convert   f 
to   kHz 

184           +                 -55 
1S5      RCLS          36  89 

1 86  -                -45 

187  R/S                51 

Compute  & 

Display  TL 

e               n 

137  X2                53 

138  STOA          35  11 

Store   f 

188     tLBLo.     21   16  11 
169       STOC          35  13 
196          x                 -35 

191  XZ'f              -41 

192  r                -24 

193  STOD          35  14 

194  QSBb     23  16  12 

195  RCLC          36  13 

196  GSBh     23  16  12 

197  -                -45 

198  x                -35 

199  RCLC          36  13 
266           t                 -24 

Subroutine 

139            .               -62 
148             i                 81 

141  x                -35 

142  1                 81 

143  RCLti           36  11 

144  +                 -55 

145  t                 -24 

1 46  1+                 56 

Attenuation 
Coefficient 

First   term 

147  4                 64 

148  8                 86 

149  RCLh           36  li 
J  58           x                 -35 
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Step  Keys   Code 


281 
282 
233 
284 
285 
286 
287 


28S 

28? 
218 
211 
212 
213 
214 
215 


RCLi 

ISZI 

hBS 


R7N 


36  45 

16"  26  46 

-24 

16  31 

82 

-35 

24 


Kxplana  tion 


A? 

Subroutine 

(Continued) 


Step   Keys    J  ode 


■L'xp  lariat  I  on 


$IBL6 

CHS 
i 

+ 

RTN 

R/S 


£1  16  12 
53 

-12 
81 

-55 
54 
24 
51 


cos   to    sin 
Subroutine 


Z   Gain   and  Transmission   Loss   Program 
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